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Abstract. Massive, large scale star counts in the J (1.25 
/im) band provided by the Deep Near Infrared Survey of 
the Southern Sky (DENIS) are used for the first time to 
draw out an extinction map of the Chamaeleon I dark 
cloud. We derived a 2' resolution map of the cloud from 
J star counts within an area of I? 5 x 3° around the cen¬ 
tre of the cloud using an adaptive grid method and ap¬ 
plying a wavelet decomposition. Possible contaminating 
young stellar objects within the cloud are removed, al¬ 
though they are shown to have a negligible effect on the 
counts. A comparison of our extinction map with the cold 
contribution of the IRAS 100 /rm emission shows an al¬ 
most perfect matching. It is shown that J star counts su¬ 
persede optical counts on Schmidt plate where Ay > 4. 
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1. Introduction 

The Chamaeleon I cloud is the most obscured region of 
the Chamaeleon dust-molecular complex. An extinction 
map of this cloud has been drawn out by Gregorio Hetem 
et al. (1988) using star counts on ESO(B) plates. It shows 
a maximum of Ay ~ 6, and the extinction profile across 
the cloud centre shows an abrupt growth between 2 and 6 
Ay , followed by a plateau at the peak of extinction which 
looks like a saturation. As stars become too scarce, opti¬ 
cal counts can no longer be used to derive the extinction 
accurately. Near infrared (1-2 /rm) star counts are more 
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appropriate to probe regions where Ay ^ 4 since an ex¬ 
tinction of 10 visual magnitudes drops to only « 3 mag¬ 
nitudes in the J band at 1.25 ^m. Massive star counts 
in the near infrared are made possible for the first time 
thanks to the large scale near infrared survey, DENIS, cur¬ 
rently in progress (Epchtein, 1997). The aim of this paper 
is to investigate in detail the extinction toward the Cham 
I cloud using this new wealth of data. 

The Cham I dark cloud is located at 6 = —16° and its 
current distance estimate is 140 pc (Whittet et al, 1987). 
Its high galactic latitude implies a small number density of 
background stars which limits the spatial resolution of the 
extinction estimation, but, on the other hand, the proba¬ 
bility of crossing several clouds on the line of sight is low. 

2. Observations 

The observations presented here have been collected as 
part of the DENIS survey between January and May 
1996 at La Silla (Chile) using the ESO 1 meter telescope 
equipped with the specially designed 3-channel camera 
(Copet, 1996). They cover an area of 1?45 x 2?94 centered 
at a = lli^Od"", 6 = -77°30“ (J2000) in three bands, 
namely I (0.8 /rm), J (1.25 /rm) and Ks (2.15 fim). They 
consist of 13 strips each involving 180 images of 12' x 12' 
taken at constant RA along an arc of 30° in declination. 
The overlap between two adjacent strips reaches 75% in 
the Cham I cloud because of the proximity of the south 
pole. Wherever a star is picked up in two adjacent strips, 
position and flux values are averaged. Limiting magni¬ 
tudes are 18, 16 and 14 at 3a in I, J and bands, 
respectively, and position accuracy is 1" in both direc¬ 
tions. 

3. Star count method 

Usually, the extinction is evaluated by comparison of star 
counts in the absorbed region and a nearby area assumed 
to be free of obscuration (Wolf diagram method). Star 







2 


L. Cambresy et al.: Extinction map of Chamaeleon I molecular cloud with DENIS star counts 


counts are performed by adding up the stars up to a 
given magnitude (or in a given magnitude range, e.g., 
± 5 ) within a grid of fixed squares. The step of the grid 
is a compromise between the stellar density and the spa¬ 
tial resolution. In other words, the spatial resolution is 
underestimated wherever the extinction is low, while in 
highly obscured areas, the content of several cells must be 
merged, in order to pick up enough stars. 

We have developed a new method to investigate the 
extinction across a cloud which consists in replacing usual 
star counts by an estimation of the local projected star 
density obtained by measuring the mean distance of the 
X nearest stars. The most important advantage of this 
method is to match the local extinction: it corresponds 
to a star count with adaptable square size. Another very 
interesting advantage of the method is to provide a map 
with white noise. Therefore, we can simply estimate the 
noise by computing the standard deviation a of the mean 
distance on a part of our map with no signal. 

We obtain a map where each point represents the 
square root of the local density. The extinction is then 
easily derived by the relation : 


Ax = - log 
a 



where a is defined by : 
log (dcp) '^ = a X mx + b 


( 1 ) 

( 2 ) 


where mx is the magnitude, dd the mean distance of the 
X nearest stars in the cloud and dcp the mean distance of 
the X nearest stars in a comparison field supposed unob¬ 
scured. We found a = 0.34. We verified that the relation 
(2) is correct up to J = 16, i.e. our limit of complete¬ 
ness. Actually, we do not need a reference field to derive 
the extinction because the available data cover 30° of lati¬ 
tude and, thus, widely exceed the cloud dimension. So, we 
plotted the density versus the galactic latitude to interpo¬ 
late the density inside the cloud for no extinction. Then 
we convert Aj into visual extinction using the extinction 
law of Cardelli et al. (1989) for Ry = = 3.1. So 

4 ^ = 0.282. 

A limitation to star counts behind molecular clouds, 
is the possible presence of young stars embedded inside 
the cloud itself. To draw out a reliable extinction esti¬ 
mate, the counts must be dominated by background stars. 
Therefore, we have attempted to remove these spurious^ 
although interesting, objects using a colour excess crite¬ 
rion. A first iteration of the extinction estimation is car¬ 
ried out without taking into account these objects. Then, 
this map is used to deredden all stars, individually. Their 
colours J — K are compared to the main sequence star 
at the cloud distance in a colour-magnitude diagram {K 
versus J — K). Then we flag each star which would have 
an extinction even greater than 5 Ay and which are likely 


4 


6 


10 


12 


14 



0 1 2 3 4 5 

(J.KIc 


Fig. 1. Colour-magnitude diagram of dereddened stars in an 
area of 4.3 deg^. Diamonds represents young stellar object can¬ 
didates (100 on 30 000 stars), the arrow corresponds to the 
extinction vector of Ay = 10 



Fig. 2. Spatial distribution of young stellar object candidates. 
Extinction isocontours at 2, 4 and 7 Ay are overplotted 


to be intrinsically very red objets. The colour-magnitude 
diagram is presented in Fig. The spatial distribution 
of these flagged stars (Fig. ^) show essentially two clus¬ 
ters near the two AO stars HD 97300 (at the north) and 
HD 97048 (in the centre). This result is in agreement with 
several studies of T-Tauri associations in the Cham I with 
IRAS (Assendorp et al., 1990) and ROSAT (Lawson et al., 
1996). We find also several red stars farther away from the 
dense cores of the cloud. They could be real objects (T- 
Tauri, red giant), or their position in the diagram might be 
unreliable because of the uncertainty on the J — K colour. 
We emphasize the fact that our criterion identifies only the 
classical T-Tauri, but not the so called weak-line T-Tauri 
which have no or little infrared excess. Nevertheless, this 
type of T-Tauri does not concentrate in the dense cores as 
the classical one does, so they have probably little effect 
on star counts. Finally, this operation removes only 0.5% 
of stars (w 100 over 30 000). Counts are, thus, strongly 
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dominated by background stars. Then, we have built an 
extinction map from the cleaned counts. 

We can consider our map as a digitized image which 
allows to use current technics of image processing such as 
the wavelet transform to restore the image and to filter 
the noise (Starck & Murtagh, 1994). We apply the a trous 
wavelet transform algorithm to split-off the image into 4 
wavelet planes. The decomposition is made by convolv¬ 
ing the image by a low-pass filtering matrix. The differ¬ 
ence between the original image and the result of the first 
convolution gives the first plane of the wavelet transform 
which corresponds to the high frequency plane. Further 
iterations of this process provide the 4 wavelet planes and 
the final smooth plane. 

Thus, we can use the high frequency plane to identify 
aberrant points and remove them in the final image in 
order to eliminate their contribution in all the planes, by 
replacing the bad pixels by the average of the surrounding 
8 pixels. We are conscious that this process might result 
in a loss of information, but less than 0.05% pixels are 
actually corrected in this way. 

Lastly, we filter each wavelet plane using the following 
method. The noise on star counts is poissonian, but tak¬ 
ing the logarithm, as defined in equation (1) changes the 
statistical properties which are no longer poissonian. A 
Poisson noise having the standard deviation estimated in 
a region of the map with no signal is simulated. Then we 
take its logarithm and we decompose this simulated noise 
into wavelet planes. The estimation of the standard devi¬ 
ations tTi on each plane allows an adaptable thresholding. 
Then we filter each plane at 

4. Results 

In order to optimize the resolution, we choose the J band 
and a number x = 20 stars to estimate the local density. 
The extinction is characterized at each point by the mean 
distance to the 20 nearest stars. In the J band the result¬ 
ing spatial resolution is 1' for low extinction [Ay ^ 1) 
and 2'30" for Ay ^ 8. In the I band we find more than 4' 
and the Kg band is not sensitive enough. The final result 
is the extinction map presented in Fig. ^ with the IRAS 
100 /im emission isocontours in which a warm contribu¬ 
tion has been subtracted (see below). This map in false 
colours results from the recombination of the 4 wavelet 
planes. 

The standard deviation of the map is ct = 1.4 magni¬ 
tude. The resolution for high extinction allows to clearly 
identify 4 distinct maxima greater than 7 which were not 
resolved by the Schmidt plate analysis. The maximum is 
Ay ~ 9.8. There is no evidence of saturation. Using the 
relation between visual extinction and Nh column density 
(Savage and Mathis, 1979) we derive the mass M of the 
cloud for Ay > 2 : 

M = agi<P A*, = 280Mq 



Fig. 3. Extinction map derived from J band stellar counts 
and cold IRAS 100 gm emission isocontours at 20, 25, 30, 35 
MJy.sr“^ (see text) 

where a is the angular size in square radian, d the distance 
to the cloud, and ^ the mean particle mass. To calculate ^ 
we consider a gas composed of 91.5% of hydrogen (H 2 -I-HI) 
and 8.5% of helium. The mass cannot be estimated in a 
straightforward way for Ay < 2 because the slow variation 
of the extinction and the lower sensitivity in the J band 
induce important consequences on the value of the solid 
angle a which delineates the absorbing region, but the 
larger uncertainty on the mass determination comes from 
the distance estimate. 

The extinction map that we have drawn out and the 
100 /rm map of IRAS have a comparable spatial resolution. 
It is therefore tempting to cross-correlate the 2 maps to 
derive some properties of the dust grains. According to 
Laureijs et al. (1991) the IRAS 100 ^m flux consists of a 
cold and a warm components which can be split off into 
two components using the 60/100 ^m colour temperature. 
Following Boulanger et al. (1997), the cold contribution 
can be written: 

l^coid(100 p,m) = 1.67 X (/i^(100 /im) — 5 x I^{60 /im)) 

Figure shows that our extinction map and the 
cold 100 /im emission are in very good agreement. 
Note that the two areas with no IRAS contour 
(11^10“, -76°34“ and 11^08“, -77°34“) are due to the 
presence of peculiar objects such as the Infrared Nebula 
(Schwartz & Henize, 1983) and very young stellar objects 
with outflows (Jones et al., 1985). This excellent correla- 
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tion suggests that the J extinction and the cold 100 /rm 
emission have the same origin, a result in agreement with 
the Desert at al. (1990) dust model which shows that the 
100 /rm emission and the near infrared extinction are both 
caused by big grains. The warm component contribute 
strongly to 100 /rm emission, but not much to the ex¬ 
tinction. A plot of the cold IRAS 100 /rm emission versus 
visual extinction is presented in Fig. 
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Fig. 4. Variation of cold IRAS 100 pm emission with visual 
extinction. Each point corresponds to a position on the map 
spaced of V in both directions 


5. Conclusion 

The extinction map of the Chamaeleon I cloud has been 
significantly improved for extinction greater than 4Ay 
with respect to previous maps obtained from star counts 
on Schmidt plates. Four distinct maxima are detected and 
we reach 9 visual magnitudes of extinction without degra¬ 
dation of the resolution. This result has been obtained 
both by exploiting the massive star counts in the J band 
provided by DENIS, and by applying a variant of the clas¬ 
sical star count method which is adapted to large varia¬ 
tions of extinction and a wavelet analysis of the extinction 
map. Moreover, DENIS give us the opportunity to inves¬ 
tigate the cloud at a larger scale than the earlier inves¬ 
tigations which were limited to small regions around the 
visible reflection nebulae. 

The comparison with the cold IRAS 100 pm is strik¬ 
ing. Each of the 3 most important extinction maxima cor¬ 
responds to a peak of cold IRAS emission. The IRAS flux 
is therefore a good indicator of extinction. We plan fur¬ 
ther investigations of the relation between the cold IRAS 
emission and extinction when new CO observations at a 
comparable spatial resolution will be available. 

Finally we stress the fact that our map has been de¬ 
rived from J star counts converted to visual extinction 
with Rv = 3.1. This value of 3.1 is the standard esti¬ 
mate for diffuse interstellar medium, but it can reach 5.5 


in dense molecular cloud cores (Whittet et al., 1987). Ac¬ 
cording to Cardelli et al. (1989) : 


A 

= 0.4008 - 
Av 


0.3679 

Rv 


So, we obtain 0.282 for Rv = 3.1 and 0.334 for Rv = 5.5. 
This means that we overestimate the extinction by a factor 
1.18 if we choose Rv = 3.1 rather than 5.5. Extinction is 
generally expressed in visual magnitudes. This choice is 
not the best because the extinction law in the visible range 
depends on the composition of the medium. On the other 
hand, the extinction law in the infrared seems universal. 
Therefore, it should be better to refer the extinction to 
near infared magnitudes. Lastly, it appears that J DENIS 
data with a detection limit of lO*** magnitude are better 
adapted to investigate the obscuration of regions where 
the extinction is larger than 4 magnitudes, rather than 
the Schmidt plate even at a usual detection limit of 22"'^ 
magnitudes. 
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